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A unified treatment for the propagation of classical waves in inhhomogeneous media is proposed. 
We deal with four kinds of waves, they are the acoustic wave in fluid, the elastic shear wave in two 
diemnsional solid, and the E- and H-polarized electromagnetic waves in two diemnsional lossless 
medium. We first show that a universal wave equation governing the wave motion of all these 
four kinds of waves can be derived. We then introduce an auxiliary field, and give the universal 
expressions for the energy densities and energy flows of these waves. This unifled treatment provides 
intuitive insights, which may be helpful in understanding the essential physics of various wave 
phenomena, or useful for designing new photonic and sonic devices. 

PACS numbers: 41.20.-q, 41.20.Jb, 43.20.+g 



I. INTRODUCTION 

Periodic media such as the photonic crystals 
0JiJlJ3 and phononic crystals [1 H Efl [H UJ, 11^ 
Il4l Il5l lla. HtL ll^ have received considerable attention 
in the past decades because of their interesting physi- 
cal properties and apphcation potential in designing var- 
ious new devices. Waves propagating inside these media 
will be modulated by the periodic structures. There- 
fore, the behavior of these waves will no longer be the 
same as they do in a free space, and the so called fre- 
quency band structures appear Besides, if we de- 
stroy the periodicity to make these structures disordered, 
then under certain conditions there might appear an- 
other interestin g ph enomenon called Anderson localiza- 
tion mmiUlMil, which is caused by the multiple 
scattering of waves in the media. Recently, the fascinat- 
ing and unusual negative refraction (NR) phenomenon 
I25I |26| in a kind of artificial media called left-handed 
materials (LHM) [23, received even more attention 
and has alr eady stimulated a l arge amount of investa- 

gations liilMlIllllSlllllEllllalllilliliS 

Even more interestingly, the newest research results have 
shown that the NR-like phenomena can also happen in 
acoustic systems s42j . 

All of these developments are about classical waves 
including the electromagnetic (EM) waves, the 
acoustic (AC) waves, and the elastic (EL) waves. Clearly 
there are similarities between the EM and AC systems 
(hereafter both the AC and EL systems will be denoted 
as the AC systems) , and in the past decade these similar- 
ities indeed stimulated some parallel developments about 
the studies on EM and AC systems, such like the photonic 
and sonic crystals. Although these similarities have been 
noticed since as early as in the 19th century ^431, however, 
up to now it seems that a systematic study and a com- 
plete understanding on them are lacking. This situation 
leads to some inconveniences. For example, there have 
been proposed useful guiding rules for designing photonic 
crystals 0, S H and phononic crystals [ll, 21, E3, HI 
with large band gaps. However, it is not easy to un- 



derstand the explicit meanings of these rules. Besides, 
since in most studies the researchers put restrictions on 
the values of the material parameters (for example, in 
a photonic crystal one usually assumes the permeability 
/i = 1), the usefulness of these rules might be very restric- 
tive. By extending the value of material parameters to 
the unexplored regions, some interesting or unexpected 
phenomena might be found 22, _23, ^4] . An abstract and 
unified treatment are thus useful and can avoid unnec- 
essary repeated works. It can also help us to find the 
essential physics that governs various wave phenomena. 
Furthermore, since this treatment unifies the AC and EM 
systems, anyone who is familiar with only one side, say, 
the AC systems, can "translate" his knowledge to the 
other side, and thus can help him to understand the ap- 
parently different systems more easily. 



In this paper, we shall show that for AC wave in 
fluid, the EL shear (ELSH) wave in two-dimensional 
solid, and the E- and H-polarized EM waves in two- 
dimensional nonabsorptive media, a unified theory can 
be constructed, which gives us a universal description of 
the dynamical behaviors of the above mentioned classical 
waves. We believe that a further development of this the- 
ory can help us much in clarifying the essential physics 
of the classical wave systems. 



This paper is organized as follows. The explicit wave 
equations for the four kinds of waves mentioned before 
will be given in the next section. A unified treatment 
for these waves is developed in section HI, in which we 
introduce an auxiliary field. Using this auxiliary field, 
not only a universal wave equation, but also all physical 
quantities such as energy densities and energy flows can 
be expressed in a unified manner. We then discuss in sec- 
tion IV the the applications of this unified theory. Section 
V is the summary of this paper. Finally in Appendix A 
we derive the conservation laws for elastic waves. 
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II. THE WAVE EQUATIONS 



Acoustic waves 



In this section we first review the wave equations of the 
EL, AC and EM waves, respectively. We then demon- 
strate that for 2D propagation, the wave equations for 
the above mentioned four kinds of classical waves have 
the same form and thus can be treated in a unified man- 
ner. 



A. Elastic waves 

Let u be the deformation of the medium at certain 
point from its equilibrium position due to the stress. In 
terms of Cartesian coordinate system, the propagation of 
EL waves in an inhomogeneous medium is governed by 



d A duj 
dxi \ dxj 



d 



dui duj 
dxj dxi 



dxj 



(1) 



where p is the mass density, Ae and are the Lame 
coefficients, and 



Tij = Ae {diui)6ij + fie [diUj -f djUi) , 



(2) 



is the ijth component of stress tensor, all quantities are 
position dependent. In terms of vector notations, the 
above equation can also be written as 



One special case of Eq. is the wave equation for AC 
waves in fluid. In fluid, we assume /^e = (no shear force) 
and thus Ct = 0. Now we denote the phase velocity of the 
wave as c, where c = c/. The pressure p and vibration 
velocity v can be derived via the relations 



du 
'dt' 



p=-AeV-U, v=— . (8) 

Substitute (O into and one flnds the wave equation 

(9) 



1 d'^p 



and the equation relating v to p: 

dw 1 
at p 



2. Elastic SH waves 



(10) 



Another special case of Eq. is the wave equation for 
ELSH waves. Suppose the medium considered is homo- 
geneous along the direction of z, and wave is propagating 
along a path lying on the xy plane, then the u vector can 
always be written as a sum of two decoupled vectors: a 
u^j, vector that lies on the xy plane, and a vector that 
parallels the z-axis. Now if u^y = 0, i.e., SH wave, then 
u = Mz and the wave equation becomes 



5t2 dx. 



In a homogeneous medium u can always be written as 
the sum of a gradient and a curl, i.e.. 



U = U; + Ut = V(f> + V X 



(4) 



where u; = V$ and Uj = V x ^ represent the longitudi- 
nal and transverse components of the waves, respectively. 
The fields $ and ^ satisfy 



1 



= 0, 



1 

V2\^ - —— = 



(5) 



P-g^ = V • (/ieVw) - V • (pc'tVu). (11) 
For the convenience of the following discussions, we define 



_ 1 _ 1 
^~ Pe pcf ' 

and thus Eq. (|11|1 can be rewritten as 
V J vet 



,2 g^2 



B. Electromagnetic waves 



(12) 



(13) 



where 



' Ae + 2pe 



(6) 



are the phase velocities of the longitudinal and trans- 
verse waves, respectively. From these results the Lame 
coefficients can be expressed as 



Now we turn to the discussion of EM waves. Adopting 
the Gaussian unit, and denoting the speed of light in 
vacuum as cq. Maxwell's equations in a source free and 
locally isotropic medium are written as 



„ ^ 1 9B pdU 

VxE = ^-ii 

Co at Co ot 



(14) 



(7) 
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and 



V • D = V • (eE) =0, V • B = V • (pU) = 0. (16) 

Here e — e{r) and /i = fi{r) are position dependent per- 
mittivity and permeability. From (|14|l and (|15|l . we can 
esaily derive the wave equations in terms of E and H 
fields: 



V X 



V X 



V X E 



V X H 



92 E 



(17) 



(18) 



1. E-polarized and H-polarized EM waves 

Like that in the discussion of the elastic waves, we 
now assume that the medium is translational invariant 
along z. It is well known that in this case if an EM 
wave propagating along a path on XY plane, then elec- 
tromagnetic fields can always be decoupled into the E- 
polarized and H-polarized waves. For E-polarized waves, 
we mean E = Ez. Remember that in this system 
9(anything)/9z = 0, thus we have 



V X 



V X E 



= -zV . 



V M / 



So, Eq. l(T7|) becomes 
V 



/V£;\ _ e d^E _ 1 d^E 
\ fi J Cq dt^ ^c^ dt^ 



Similarly, for H-polarized waves, we have 



c2 dt^ 



H 



(19) 



(20) 



(21) 



Here c = co/y^ represents the speed of light in the 
medium. 



will show in the later discussion that the complete unified 
description can be constructed using the wave equation 
of the auxilliary field. 

Consider an auxiliary field satisfying the wave equation 



a 



1 



0, 



(23) 



which has the same form as in Eq. 1)22(1 with a = a(r) and 
c = c(r) being position-dependent parameters. Through- 
out the paper we assume d^jdz = 0, i.e., $ = ^[x,y,t). 
Under this assumption the system becomes effectly two 
dimensional. We define the dynamical variables as ip and 
Q: 



Q 



-v$. 



(24) 



We will see in the later discussion that $ plays the role 
of the vector potential in the electromagnetics. Using the 
difinition (|24ll the wave equation l|23|) can be replaced by 



5Q 

'dt 



1. 



dip 
'dt 



-ad-V ■ Q. 



(25) 



Note that if we choose the Lagrangian density 
2 



1 /9£ 
2^ Vdt 



1 



2ac2 



-Q2, (26) 



then Eq. (|23|l and (|25|l can be deduced from the Euler 
equation of motiom 



9$ ^ a5 



9(V$) 



0. 



(27) 



Here we have employed the notation $ = d^/dt. With 
this lagrangian, the canonical momentum H and the 
Hamiltonian 7i are given by 



H 



9$ 



1 



1 



(28) 



and 



III. UNIFIED TREATMENT FOR CLASSICAL 
WAVES IN 2D SYSTEMS 

All the wave equations of the four cases discussed be- 
fore have the same form 



\ a 



1 d'^U 



(22) 



thus their wave equations can be treated in a unified man- 
ner. However, this is not enough if we want to construct 
a complete unified description. Such a theory must also 
deal with the quantities like energy density and energy 
flow in a unified manner. To achieve this goal, we now 
introduce a field whcih we call the auxilliary field. We 



n 



H$ - £ 



-H^ 



a 



2ac 



1 

2^ 
1 



2 



(29) 



A. The Energy conservation law 



In terms of Eq. (|23|) , we now derive the energy conser- 
vation law. Multiplying — $ on Eq. (|23|l . we have 



/ ^v* \ v<i> • v$ $ a$ 

lay a at? dt 
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a J dt 



1 



1 



2ac2 ^ ' 2a 
Define the energy flux J and energy density lA as 

a 



= 0. (30) 



1 



1 



and 

tlien Eq. (|30|l becomes 



^,2 



a 



(31) 



(32) 



(33) 



which is the energy conservation law (the "continuity 
equation" for the energy). Comparing Eq. H32|l with 129|l 
confirms that lA is indeed the energy density. 



B. Identification of J and U with Icnown physical 
quantities 

1. J and U for acoustic waves in fluid 

For the acoustic wave propagation in 2D we define 



$ = / pdt, 0L = p, 



(34) 



then 



<^ = P, Q = v. (35) 
The energy flux J and energy density U are given by 

J=pv, = ipi,2 + _l_p2^ (36) 

which are consistent with the results from the 3D formula 
(See appendix A). 



2,. J and U for elastic SH waves in solid 



For the SH waves u — u{x,y,t)z. We define 



1 1 
<P = u, a = — = — 5-, 

Me pet 



(37) 



then 



ip = u = V, Q = —fipVu = — T • z, (38) 



where v = represents the vibration velocity of the 
media and T stands for the stress tensor. 

Now we denote n (||Vm) as the direction of Q, i.e., 



Q = Qn, Q=|Q| 



(39) 



Suppose there is a region enclosed by a surface with nor- 
mal vector n (here we have n±z). The surface exerts 
a shear force fs (along z.) on its surrounding medium, 
which is given by 



fs = -Tn = Qz. 



(40) 



From this observation we find that Q in this case is a 
vector that has the magnetudc of the shear force and the 
direction of Vw. The energy fiux J and energy density U 
are given by 

J^-^ieWu, U^^v^ + ^{Vuf, (41) 

which are consistent with the results from the 3D formula 
(See appendix A). 



3. J and lA for E-polarized wave 

For EM waves, we start with the consideration of the 
E-polarized wave. In this case the dynamical variables 
are E and H fields. Since V • B = 0, in terms of the 
vector potential A we have B = V x A. Now we choose 



A = -VIttcq $z 



and 



then we have 



E = 



1 OA 

Co dt 



co/e, 



47r 
Co 



H = -V X A = 



Q X z. 



(42) 
(43) 

(44) 
(45) 



Using these relations then Eq. H25(l indeed gives the cor- 
rect dynatimcal equations for E-polarized waves (i.e., 
Eq. (|14|) and (|15|) '). In addition, the textbook formulas 
for the Poynting vector 

S = -^E X H = (^z X (Q X z) = <^Q = J, (46) 
47r 



and the energy density 



2ac2 



^p^=U. (47) 



give exactly the same results as those obtained from our 
formulas for energy flux (Eq. (|31|) ) and for energy density 
(Eq. 1S3). 



4- J and U for H-polarized waves 

In this case the dynamical variables are also E and 
H fields. Since we consider only the source free case, we 
have V • D = 0. Thus we can introduce a vector Am such 



5 



that D = V X Am- Here A™ represents another kind of 
vector potential, used in deriving the D field. Now we 
relate to $ via 



and choose 

then we have 
H 



a = e/co, 



1 dAr, 
Co dt 



ac 



co//i, 



Co 



1 . /47r _ , 

E = -V X A„i = —\ — Q X z. 
e V Co 



(48) 
(49) 

(50) 
(51) 



Using these relations then Eq. H25I) gives the correct dy- 
namical equations for H-polarized waves (i.e., Eq. 114|) 
and l|15(l ). Furthermore, the Poynting vector 



co_ 
47r 



H X E = —ipz X (— Q X z) ~ tpQ, 



(52) 



together with the electromagnetic energy density 



ip^^U. (53) 



are the same as those obtained from Eq. H31|) and 
Eq. (EH). 



The time averaged energy flux and energy density are 
defined as 



1 f 
^ -j^ 3{t)dt 

= iRc(^*Q 
2ac \ ik 



(59) 



and 



1 
T 



U{t) dt 



1 / 



4ac' 



(60) 



Monochromatic waves 



Here T = 2tt/uj is the time period, and k = A:(r) = Ci;/c(r) 
is the wave number in the media. 



IV. APPLICATIONS 

Now we turn to the applications of this unified theory. 
In the following subsections we will discuss the use of the 
theory, only briefly, on two topics: the wave scattering in 
a two dimensional system, and the band gap engineering 
problem. 



Now we consider the monochromatic wave. The aux- 
iliary field $ can be written as 




Here $ cx e~*"* satisfies <& = —iuj^. is the complex 
conjugate of 

The ip and Q now become 

>f = \{p=+r), q = ^(q + q*), (55) 

with (p and Q given by 

(p^-iuj^, Q = --V$. (56) 
a 

Since both (p and Q contain a factor e^*'^*, they can be 
factorized as 

(^(r,t)=V(r)e-^"*, Q(r, i) = f *, (57) 

where ip is a, time independent field. 

The wave equation Eq. H23() now becomes 

\ a(r) / a(r)c^(r) 



A. Scattering in a two dimensional system 

We discuss only a simple case. Suppose a (monochro- 
matic) incident wave from a source propagates to and be 
scattered by a circular cylinder located at the origin. 

Now we assume the radius of the cylinder is a. The 
material parameter a{v) inside and outside of the cylin- 
der are given by constants ai and a, respectively. Sim- 
ilarly, the wave speed are given by ci and c. We also 
define fci = w/ci and k = to / c as the corresponding wave 
numbers, and denote the wave inside and outside of the 
cylinder as ^/'i(r) and ■0(r). According to Ea. (|58|l . they 
satisfy the Helmholtz equations, 

(V^ fci)!/-! = 0, (V^ + fc^)7/' = 0. (61) 

Since the Bessel and Hankcl functions of all orders 
{J„, Hn\n = — cxD — > cxd} forms a complete set of eigen- 
functions, the solutions of ^ and tpi can be written as 
the sums of them: 

oo 

n— — oo 

oo 

Mr,0) = J2 + A.Jn(fcir)]e™«(62) 
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To determine the coefficients, we have to find the 
boundary conditions. According to Ea. H58|) . and by using 
the divergence theorem, we get the boundary conditions 
at r = a: 



Here 



1 aV'i 



ai dr 



1 

a dr 



(63) 



Furthermore, since the wave source is outside of the 
cyhnder, the wave ampHtude inside the cyhnder cannot 
goes to infinity, we therefore have C„ — 0. Substitute 
Eq. (IH2Jl and C„ = into , we have 

A„H^P{ka) + BnJnika) = DnMha), 

-[A„Hi'^' (ka) + BMka)] = — A.4(fcia). (64) 
a ai 

These two equations lead to 

An ghJ'n{ka)Jn{ka/h) — Jn{ka)J'n{ka/h) 



Bu H^n\ka)J'n{ka/h) - ghni^^' (ka)Jn{ka/h) 
and 



Dn 



2gh/ (inka) 



Bn H^}\ka)J'n{ka/h) - ghH^n^' {ka)Jn{ka/h) 



(65) 



(66) 



Here we have defined g = ai/ a and h = ci/c. In deriving 
Ea. (|^ . we have used the identity 



2i 

TTX 



J„(a;)ffW'(x)-4(a;)ffW(x) 
Suppose the incident wave is a plane wave: 

oo 

V'o(r) = e"'"'" = e'fc'-™«» = ^ ^"J„(fcr■)e™^ (67) 

71— — OO 

then we have 

Bn = i". (68) 

Now the coefficients An and Dn can be obtained from 
Ea. (|^ and (|^ . The scattered wave is then given by 



V'.ca*(r) = V(r)-V'o(r)- ^nH^^Hkr 



(69) 



If the incident wave is not a plane wave, the only differ- 
ence will be to give a different set of -B„'s, and calcualte 
its corresponding An and -D„'s from Ea. (|65() and (|66|l . 



B. Band gap engineering 

Now we turn to the discussion about the photonic and 
phononic crystals. In these wave crystals a(r) and c(r) 
are periodic functions of r, i.e.. 



a(r + R)=a(r), c(r) = c(r + R) 



(70) 



R ~ niai + n2a2 



(71) 



is an arbitrary translation vector, ni and n2 are two in- 
tegers, and ai and sl2 are the fundamental translation 
vectors 0|. According to Block's Theorem, the eigen- 
wave in wave crystal has the form 



V'K(r)=e'^-^(r) 



(72) 



where K is the Block wave vector, and ^(r) is a periodic 
function, satisfying 



^(r)=^(r + R) 



(73) 



The band structure of the wave crystal can be cal- 
culated in a unified manner from Ea. (|58|l by using the 
method that used in the second paper of Ref. We 
will not go through the details of the band structure cal- 
culations but instead give a brief discussion about the 
band gap engineering. 

The main aim of the band gap engineering in the wave 
crystal study is to create a large band gap. Recently, 
Zhang et.al have developed a systematic method for en- 
larging a photonic or phononic '17*1 band gap. Their 
methods are based on the variational principle. Their 
formulas tell them how and where to vary the material 
parameters can enlarge a band gap that already exists. 
Here we will derive the unified version of their formulas 
(See Eq.(5) and Eq.(7) of Ref? and Eq.(2) of ReflH. 

Suppose V'K(r) = V'(r) is an eignfunction of the oper- 
ator 



-a(r)c^(r)V 



1 



a{r) 



and = Lu'^ is the corresponding eigenvalue. That is, 
tp{r) satisties Eq. Multiplying Eq.jSSj) by ■0*, we 

find 

V.r^V^ + ^ = 0. (74) 
\ a J a ac^ 

According to Ea. (|72|l . the quantity 

a a 

is a periodic vector field. Integrating Ea. H74() over a unit 
cell, we get 



dV^^^ / /3(r)|V(r)|2dV. (75) 

cell Cy.\J^) Jcell 



Here we have defined 



Q;(r)c^(r) 



Rewriting Ea. l(58|) as 



(76) 



(77) 
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we can derive 



V • 



5 ( ^ ) VV' 



(78) 



Note that in any case we do not vary the Bloch wave 
vector K. 

Muhiplying Ea. lfTI^ by i/'* and integrating the result- 
ing product over one unit ceU, we find 



cell 



a J a 

2 j2„ 



Si 



cell 



(79) 



Multiplying the complex conjugation of Ea. (|77(l by 5ip, 
we can also derive 



Pij^*5ijj.d^r 



(80) 



cell 



By substituting Eq.||Sni) into Eq.lO, and using Ea. ((75|) . 
we finally obtain 



5lj 

LU 



Liim'd'r 



(81) 



This is the unified version of the band gap engineering 
formula derived by Zhang et. al. 0, IT^ . 

For acoustic wave we have if) = a = p and /? — 
1/ac?, and hence 



5p 



pj (pc2)2' 

We therefore have from Eq.JHIJ and Ea. l75|l the result 

L 



b(r)|-'^(p(r)c^(r)) ,2 
Icell (p(r)c2(r))2 " ' 



|p(r 



^cell /9(r)c2(r) 
^cell P^{^) 



dh 



' f 



cell p(r)c^(r) 



d2r 



(82) 



which is just the Eq.(2) of Ref 17. Resuhs for the EM 
waves can also be obtained in a similar way. 

We now give several remarks about Eq.jTSJ. 

A. In fact, Eq. (|75|l is a consequence of the harmonic 
property of the wave field. In the second line of Ea. l6()|l 
we note that the energy density can be written as the sum 
of two terms: a|Qp/4 = |VV'p/4atj2 and |v5p/4ac^ = 
\tp\^ /Aac^. Let's call them the type I energy and the type 
II energy. Eq.(J7SJ) means that the time average of this 
two kinds of energies in one unit cell is equal. This is 
reasonable because to form a harmonic wave the type I 



energy has to transform itself to type II completely and 
then transform back in one time period T/uj. 

B. From Ea.lfT^. we have 



Jcell " '"^ 



a(r) 



(83) 



which yields 



Inc^^ 



In 
In 



ai a{r) 
I3{r)\ij{r)\^d'r. 



(84) 



cell 



Taking the variation of Ea. (|84|l and keeping the \tlj\ and 
jVi/'l terms unchanged then gives us Eq.JHIJ. Why can 
we assume in this variational procedure that \ip\ is un- 
changed? The reason is that the change of IV'I caused by 
the first order perturbation of 1/a and (3 is second order. 

C. Although in a unit cell the total type I energy is 
equal to the total type II energy, however, their distri- 
butions are different. Consider a binary system, i.e., the 
wave crystal that composed of two kinds of homogeneous 
materials. A unit cell of such a wave crystal can be de- 
vided into region a and region 6, with material parame- 
ters Oiajda and Of,, /3b, respectively. According to Ref. 
IE IE Q> it is the impedance (3 /a that determine the 
size of the band gap. More explicitly, a wave crystal 
with large (iajo-a and small Ph/ah-, or a small Pa/oia and 
large /3b/ab, will most probably have a large band gap. 
This seems to imply that to have a large gap we have to 
design the a and /3 parameters such that the type I and 
type II energies are efficiently separated in the unit cell. 



V. SUMMARY 

In this paper, we have shown that for four kinds of 
classical waves propagating in two-dimensional media, a 
unified treatment can be constructed. These waves are 
the AC wave in fiuid, the EL shear (ELSH) wave in solid, 
and the E- and H-polarized EM waves in nonabsorptive 
media. This unified theory helps us to find the essen- 
tials of various wave phenomena and thus give us a more 
complete and comprehensive understanding of these phe- 
nomena. 
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APPENDIX A: CONSERVATION LAWS FOR 
ELASTIC WAVES 

The equation of motion for elastic wave is of the form 
pil^^djTj^. (Al) 
Multiplying both side by Ui, we find 

dt ( puV2 ) = d, (T,, iii) - T,, dt {djUi), (A2) 

where 

Tji dt {djUi) 
= [ Ae(V • u)(5y + fieidiUj + djUi) ] dt (djUi) 
= Ae (V • u)at (V • U) 

+ (^e/2) {djUi + diUj)dt{djUi + diUj) 
= A, (V • u)dt (V • u) + E.jdtEj, 
= dt [Ae (V • uf/2 + ^xe Tr (E^)] . (A3) 

In deriving l|A3|l we have defined a symmetric matrix E: 



(A4) 



Substituting (|A3(I into (|A2|I one finds the conservation 
law 



For the special case that acoustic waves propagating 
in fluid, we have /ie = 0. Using the relation p = — AV • u, 
we find 



(A7) 



where I is the unit matrix. 

Denoting u as v and remembering that Ae — pc^ then 
one finds 



J=pv, Z^ = |v^ + _p2. (A8) 



Another special case is the elastic SH wave. In this case 
u — uz. Remember that du/dz = 0, therefore V • u = 
and we have 



T 
E 

Tr(E2) 



2/ieE 

1 



(A9) 







^ \ OU OU Q 

dx dy 



(AlO) 



V • (T • li) 



d_ 

dt 



^u2 + ^(V.u)2+MeTr(E2) 



Thus the energy flux J and energy density U are 



(A5) 



J^-T-ii, Z^ = ^u2 + ^(V.u)2+/,,Tr(E2). (A6) 



The energy flow and energy density are thus given by 



J = -/iet'V«, ^v2 + ^(Vm)2. (All) 



Here v = u = -u z is the vibration vecocity of the media. 
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